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We present a theory of mesoscopic fluctuations of g tensors and avoided crossing energies in a 
small metal grain. The model, based on random matrix theory, contains both the orbital and spin 
contributions to the g tensor. The two contributions can be experimentally separated for weak spin- 
orbit coupling while they merge in the strong coupling limit. For intermediate coupling, substantial 
correlations are found between g factors of neighboring levels. 
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I. INTRODUCTION 

Recent developments in nanofabrication techniques 
have allowed for the resolution of individual "particle- 
in-a-box" energy levels in small metal grains or semicon- 
ductor quantum dots using tunneling spectroscopy.Efcl In 
the absence of a magnetic field, the energy levels e M are 
two- fold degenerate (Kramers' degeneracy). An applied 
magnetic field B lifts the degeneracy; the splitting of the 
doublet is described with the help of a u g factor" , 

Se» = VbqB, (1) 

where /is = is the Bohr magneton. A cartoon of the 
magnetic- field dependence of the energy levels is shown in 
Fig. [lj Whereas g = 2 for electrons in vacuum, in a metal 
grain the g factor can be different from two as a result of 
spin-orbit scattering. Recently, the magnetic-field depen- 
dence of particle-in-a-box leve Is in metal grains have been 
measured by two groups Su Measured g factors range 
from 0.1 to 2, depending on grain size, material, and, in 
the case of Ref. [|, doping with heavy ions. 

Unlike in bulk metals, where g factors are used to de- 
scribe the effect of spin-orbit coupling on the band struc-. 
ture, g factors in a metal grain are not a "bulk" property.B 
Not only does the typical value of the g factors depend on 
the size of the metal grain, g factors also depend on the 
microscopic details such as the impurity configuration, 
the location of defects, and the form of the grain bound- 
ary. As a result, different energy levels in a metal grain 
have different g factors. Moreover, even if the metal grain 
is roughly spherical and without lattice anisotropy, the 
presence of impurities breaks the rotational symmetry 
on the microscopic scale, causing g factors to depend on 
the direction of the applied magnetic field. A statistical 
description of the level-to-level fluctuations of g factors 
in metal grains has been formulated by Matveev et al.hQ 
and by Halperin and two of the authorali using random 
matrix theory (RMT). Petta and RalphEl measured g fac- 
tors for up to 9 consecutive levels in nanometer-size Cu, 
Ag, and Au grains and found good agreement with the 
distributions of Refs. |lo|,^T]. The dependence on the di- 
rection B of the magnetic field is taken into account by 
replacing the g factor by a " g tensor" Q,\3 




Magnetic Field (B) 

FIG. 1. A cartoon showing the definitions of the g-factors 
and the avoided crossing energy A. At zero magnetic field, 
all energy levels are doubly degenerate. A magnetic field 
splits these doublets. The ^-factor measures the size of the 
splitting of a doublet e M as a function of magnetic field, see Eq. 
(III). The avoided crossing energy A is the minimum distance 
at the first avoided crossing of neighboring energy levels, see 
Sec. [Fv|. 

fe M = /i B B(5 T ^B) 1/2 . (2) 

(The g tensor carries a subscript /i to reflect its depen- 
dence on the energy level e M .) The g-factor ([j]) for a 
magnetic field in the z direction is the square root of 
the tensor element Q zz . A measurement of, full g tensors 
in Cu grains was reported quite recently.cl Again, good 
agreement was found between the experimentally mea- 
sured <?-tensor distribution and RMT. 

The effect of the spin-orbit interaction on the wave- 
functions in a metal grain can be described by a dimen- 
sionless parameter A, 

T so 

where r so is the spin-orbit scattering time and 8 is the 
mean spacing between Kramers' doublets in the grain 
(in the absence of the magnetic field). The effects of 
spin-orbit scattering are weak if A <C 1. In that case, 
wavefunctions are real and have a well-defined spin; the 
electron magnetic moment is close to its vacuum value 
g = 2. In the opposite limit of strong spin-orbit scat- 
tering, A > 1, wavefunctions are complex and have no 
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well-defined spin. Hence, the spin contribution to the 
electron's magnetic moment is strongly suppressed, com- 
pared to the case of electrons in vacuum. However, in 
addition to a contribution from the electron's spin, there 
may be a significant orbital contribution to the magnetic 
moment carried by a single electron if spin-orbit scat- 
tering is present ij-jyavefunctions are complex, and hence 
current-carrying .Ej 

Experimental estimates of A are close to zero in Al and 
range from 0.7 in a small Cu grain (S fii 0.7 meV) to 13 
in a larger Au grain (S 0.1 meV).EI A full theory of 
the combined orbital and spin contributions to the g ten- 
sor was-developed for the asymptotes A -C 1 and A ^> 1 
only.L30 Both theories calculate distributions normal- 
ized to the average {{g 2 )) 1 ^ 2 ■ In addition, Matveev et al. 
calculate both spin and orbital contributions to ((g 2 )) 1 ^ 2 , 
while Ref. [ll] considered the spin contribution only. The 
case of intermediate A, necessary for a quantitative com- 
parison with the experiments of Ref. |7| was studied in 
Ref. [ll] using numerical diagonalization of a random ma- 
trix model with variable spin-orbit scattering strength, 
but without inclusion of the orbital contribution to the 
magnetization. 

In this paper we construct a random matrix theory 
that describes both spin and orbital contributions to the 
electron g tensor. In the limit A » 1 our model re- 
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produces the g tensor distribution found in Refs. ^ 
but it also provides a simple model to numerically ob- 
tain the full g tensor distribution for arbitrary spin-orbit 
scattering strength. In addition to the distribution of 
the g tensor we also look at the correlator of g tensors 
of neighboring levels. While g tensors are not correlated 
for A = and, as we show here, for A > 1; we find 
that correlations can be substantial for A of order unity. 
The random-matrix model is formulated in Sec. |fi|; the g 
tensor distributions are considered in Sec. |l|. 

In addition to the g factors, which describe the 
magnetic-field dependence of the energy levels at very 
small magnetic fields, Salinas et al. obtained additional 
information on the magnitudes of spin-orbit scattering 
matrix elements from avoided crossings of energy lev- 
els at higher magnetic fields: For weak spin-orbit scat- 
tering, the minimal energy separation A in an avoided 
crossing between the downward moving level £ M +i,- and 
the upward moving level £ M .+ is twice the matrix ele- 
ment of the spioi-orbit coupling between the correspond- 
ing eigenstateSjB see Fig. |l| 



IV we calculate 



Fig. yj. In Sec. 
the avoided crossing energy A from the random matrix 
model, and find its statistical distribution and depen- 
dence on the direction of the magnetic field B. 



II. RANDOM MATRIX MODEL 

In this section we formulate a random-matrix model 
that describes the magnetic-field dependence of energy 
levels in a metal grain with spin-orbit scattering, taking 



into account both the Zeeman and the orbital effects of 
the magnetic field. Following the basic premises of ran- 
dom matrix theory, we replace the Hamiltonian of the 
metal grain by a 2N x 2N matrix H, 



H(X) — Hqoe 



A 



H, 



GSE 



H 



B- 



(4) 



The first two terms on the right hand side of Eq. (Q) 
describe the Hamiltonian in the absence of the mag- 
netic field; the last term Hb describes the effect of the 
magnetic field. Without the magnetic field, Tt is taken 
from an ensemble that interpolates between the Gaus- 
sian Orthogonal and Gaussian Symplectic ensembles of 
random matrix theory. The Gaussian Orthogonal En- 
semble (GOE), which is relevant for metal grains without 
spin-orbit scattering, consists of real symmetric N x N 
matrices with independently and Gaussian distributed el- 
ements, multiplied by the 2x2 unit matrix I2 in spin 
space, 



H GOE = S® 1 2 , P(S)oce 



-(TT 2 /4NS 2 )trS T S 



(5) 



Here d is the mean level spacing in the metal grain (i.e., 
the mean spacing of the Kramers' doublets). The Gaus- 
sian Symplectic Ensemble (GSE), which describes metal 
grains with strong spin-orbit scattering, consists of self- 
dual quaternion matrices.EfJ A Hamiltonian taken from 
the GSE can be parameterized as 



H GSE = i ^A <g> 1 2 + i ^2 A i ® a 3 ■ j 



(6) 



where Aq is a real symmetric N x N matrix and the Aj , 
j = 1,2,3, are real and antisymmetric N x N matrices. 
The four matrices Aq, Ai, A2, and A3 have independently 
and Gaussian distributed elements, 



-(77- /4iV<5 )tr A Aj 



, j = 0,1,2,3. 



(7) 



The crossover parameter A describes the strength of the 
spin-orbit scattering in the Hamiltonian of Eq. ([!]). The 
cases A = and A — ► 00 correspond to the GOE and 
GSE, respectively. 

The effect of the magnetic field B = (B\, B2, B3) is 
described by the term Hb in Eq. (El), 



3=1 



(8) 



where the 2N x 2N matrices Mj (j = 1,2,3) are given 
by 



" 1 sVn 



(9) 



where the Xj, j — 1,2,3, arc real antisymmetric matri- 
ces, with independent and Gaussian distributions, 
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-(TT 2 /4.N6 2 )trXjXj 



(10) 



The first term in Eq. (g) describes the coupling of the 
magnetic field to the electron spin; the second term, 
which is diagonal in spin space, describes the coupling 
of the magnetic field to the orbital angular momentum. 
The second term in Eq. (^|) was originally proposed by 
Pandey and Mehta to describe the orbital effect of a time- 
reversal symmetry breaking magnetic field on the statis- 
tics of energy levels. Ej£3 For a diffusive spherical grain 
with radius R, mean free path I, and-ieffective electron 
mass m*, the coefficient rj is given bycH 



2 / / *\2 ' 

V = {m/m ) — , 



(11) 



whereas for a ballistic sphere with diffuse boundary scat- 
tering, one haaia 



rj 2 = (m/m*) 2 — . 



(12) 



At the end of the calculation, the limit N — > oo is taken. 
Without the orbital term, the Hamiltonian Tt of Eq. (0) 
is the same as the random-matrix Hamiltonian used by 
Halperin and two of the authors in Ref. [ll]. 

The g tensor Q and the avoided crossing energy A 
will be expressed in terms of matrix elements involv- 
ing the eigenvectors of the Hamiltonian Eigenvec- 
tors ipp of the Hamiltonian (Q) are 2N component com- 
plex vectors. Their elements are denoted as ip^ (n, a) , 
where n = 1,...,N refers to the "orbital" degrees of 
freedom, and a = ±1 to spin. At zero magnetic field, 
all eigenvalues of the Hamiltonian (||) are twofold degen- 
erate (Kramers' degeneracy): each eigenvalue (fi = 
1, . . . , N) has two orthogonal eigenvectors Vv an( l 
where Ttp(n, a) = aip*(n, —a), is the time-reversed of ip. 
In the GOE (A = 0, B = 0), the eigenvectors -i/v an d 



Ttpp can be chosen such that +1) = — ^"Vv( n > ~ 1) 

is a real number and Vv( n ' —1) = T^ Al (n, 1) = 0. In that 
case, the nonzero elements ^(n, +1) are independently 
and Gaussian distributed with zero mean and with vari- 
ance i/nB (Of course, any linear combination of ^ and 
Tip^ forms a valid pair of eigenvectors for the eigenvalue 
e M as well.) In the GSE (A — > oo, B = 0), the elements of 
ipp are complex numbers with independent and Gaussian 
distributions with variance 1/2N. In both the GSE and 
the GOE different eigenvectors are statistically uncorre- 
cted. 

In the crossover between GOE and GSE, the eigenvec- 
tor distribution is more complicated than in each of the 
two basic ensembles. Unlike for the cases of the pure 
GOE and GSE, eigenvectors at different energy levels 
are correlated, so that it is no longer sufficient to look at 
the distribution of one eigenvector alone.Ej Since orthog- 
onal invariance is preserved throughout the GOE-GSE 
crossover, the problem of finding the (joint) distribution 
of one or more eigenvectors in the crossover ensemble can 
be simplified by considering their orthogonal invariants 



first. For each pair of eigenstates ?/v an( l V^j the invari- 
ants are four quaternion numbers p° v , j = 0,1,2,3. If 
we diagonalize H, writing 



H(B = 0) = U(E®t 2 )U\ 



(13) 



where U is the symplectic eigenvector matrix and the 
N x N diagonal matrix E contains the eigenvalues e M on 
the diagonal, they are 



P% 



(pjp+- 



3 = 1,2,3. 



The p 



>^j, satisfy a criterion of anti-hermiticity, 



pU 



-UJ, j — 1, 2, 3. 



(14) 
(15) 

(16) 



The orthogonal invariants p ® express orthonormality of 
the eigenvectors ip^ and Ttp^ . The remaining orthogonal 
invariants p 3 ^ are characteristic for the crossover and de- 
termine to what extend spin-rotation symmetry has been 
broken. In the GOE, we have: J2k t r (/ Ji/i cr fc)t r (p^ M cr fc) — 
4<5jj, while p£„ = if p ^ v; in the GSE, pi = for 
all p and v. An average involving different eigenvectors 
is then calculated in two steps: First, eigenvector ele- 
ments have a Gaussian distribution with zero mean aud, 
with variance determined by the orthogonal invariants. 113 
In spinor notation, where ip(n) denotes the 2-component 
spinor with elements ip(n, +1) and ip(n,— 1), these vari- 
ances are 

(VvH f Vv(™-)) = -Tf-'W) 
«(VvW t(J jVvM) = ^j 1 (pi»)++, 

(Vv( n ) Tf7 2Vv(n)) = 0, 



(17) 



With the help of Eq. ( |l7| ) any average over eigenvectors 
can be expressed in terms of the orthogonal invariants 
involved in the problem. 

What remains is to find the average over a small num- 
ber of orthogonal invariants. To our knowledge, a so- 
lution of this problem exists for the limits A <C 1 and 
A ^> 1 only. For strong spin-orbit scattering, A ^> 1, the 
solution takes the form of a surmise equating the distri- 
bution of the o-?,,. for the 2N x 2N crossover Hamiltonian 
(m to the distribution of the same-quantities for a GSE 
Hamiltonian of a smaller size 2N'j£3 

N' = X 2 N(\ 2 + 2N)/(\ 2 + N) 2 

->2\ 2 ifjV-KX), (18) 
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provided the energy difference |e M — e v \ <C X 2 5. This 
means that the elements of the matrix p> are uncorrelated 
and that they have a Gaussian distribution with variance 



i 

2N' 
1 + S 



[US 



2N 



(19) 



A similar surmise was proposed in Ref. |l5| for the eigen- 
vector statistics in the crossover between the GOE and 
the Gaussian Unitary Ensemble of random-matrix the- 
ory. The motivation of this surmise becomes clear once 
we consider the crosspver Hamiltonian (Q) in the eigen- 
vector basis of TJqoe-" In this basis, eigenvectors of the 
crossover Hamiltonian are "localized" : they are mainly 
built up from eigenvectors of Hqoe with energies inside a 
window of size ~ N'6 (with N' to be determined later). 
Since changing to the GOE basis does not change or- 
thogonal invariants, we can calculate the /9^„ using an 
effective 2N' x 2N' Hamiltonian that contains the 2N' 
relevant GOE eigenvectors only, if |/x — v\ <C N'. As the 
spin-rotational symmetry breaking term is large for the 
effective Hamiltonian, its distribution is that of the GSE, 
not a crossover. The exact relation ([l8]) between TV' and 
N is found matching the distributions of a single orthog- 
onal invariant cP„„ in the crossover Hamiltonian and in 

In the following two sections, the random matrix model 
(||) will serve as a starting point for analytical calcula- 
tions of the g tensor distribution and avoided crossing en- 
ergies in the regimes of weak spin-orbit scattering, A -C 1 , 
and of strong spin orbit scattering, A > 1, and for nu- 
merical calculations of the g-tensor distribution in the 
crossover regime A ~ 1. The case of weak spin-orbit 
scattering can be treated using perturbation theory in A; 
for strong spin-orbit scattering, we use the full eigenvec- 
tor distribution of the GOE-GSE crossover Hamiltonian 
and the surmise for the orthogonal invariants that was 
discussed in this section. 



III. STATISTICS OF THE G TENSOR 

A typical plot of the magnetic field dependence of en- 
ergy levels is shown in Fig. pi A magnetic field B = BB 
splits the Kramers' doublets e M into pairs e^.± that de- 
pend linearly on the magnitude B of the magnetic field, 



= £„ ± -de 



(20) 



with Se^ expressed in terms of the g tensor as in Eq. 
(H) above. Following Ref. 11 , the g tensor can be written 
as 



Q = G T G, 
where the 3x3 matrix G has elements 



(21a) 



Ms 

G 2j = — Im (VvIMjITVv.) , (21b) 
Ms 

G 3j = —{%\M } \%), 
Ms 

where Mj is defined in Eq. (^]), ip^ is an eigenvector of 
H at B = with eigenvalue and Ttp^ is its time- 
reversed. 

The tensor Q has three eigenvectors and three eigen- 
values <? 2 , j — 1,2,3. The eigenvectors are referred to as 
"principal axes" , the eigenvalues g\ , gi , and 173 as "princi- 
pal ^-factors" . The three principal g factors describe the 
splittings of the doublet for magnetic fields along each of 
the three principal axes. We describe the distribution of 
the g tensor in terms of the distributions of its eigenvec- 
tors (the principal axes) and eigenvalues (the principal 
g- factors). For a roughly spherical grain, the principal 
axes will be oriented randomly in space. Hence, it re- 
mains to find the distribution of the three principal g 
factors g^,2, and 3^,3. We will now consider the 

cases of weak and strong spin-orbit scattering separately. 



A. Weak spin-orbit scattering 



To leading order in A -C 1, TjX <C 1 the g tensor reads 



v yfiV 



4 Si. 



^—^ p.. — p.. 



NS 



^3 



(22) 



Here e M and e v are eigenvalues of the Hamiltonian (Q) 
at zero magnetic field and without spin-orbit scattering, 
and Aj V and Xj V are the matrix element of the matrices 
Aj and Xj between the corresponding eigenvectors 
and \t/) v ) oiH, respectively, cf. Eqs. (g) and @. In Eq. 
( p2| ) we neglected terms of order A 2 /;, (A77) 2 , which are 
small compared to Xrj, and of third (and higher) order in 
A. The second term in Eq. ( p2fl corresponds to orbital 
paramagnetism, and is of order r)X because the orbital 
contribution couples to complex wavefunctions that oc- 
cur for A ^ 0. The A 2 term is a reduction of the Pauli 
paramagnetism caused by interaction with other energy 
levelS|-and for the case of i = j = 3, agrees with earlier 
work.Ej 

The distribution of Q without the orbital contribution 
(second term in Eq. (22)) was studied in Refs. |io| and [Tl|. 
For very small spin-orbit scattering, however, the orbital 
part is found to dominate the g tensor fluctuations, since 
it is of first order in the spin-orbit scattering strength A, 
while the Zeeman contribution is of order A 2 . Whereas 
the Zeeman contribution always gives g factors smaller 
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than two — the last term in Eq. (|22|) is negative defi- 
nite — the orbital contribution can be of arbitrary sign, 
allowing for principal g factors larger than two. To il- 
lustrate this feature, we calculate the tails of the joint 
distribution P(gi, 32, 33) of the three principal g factors. 
The distribution of the tails is dominated by events where 
the spacing between the level and one of its neighbors 
or e^-i is exceptionally small, of order A<5 or Xi]5 
(whichever is larger). Hence, the tails of P{g ll g 2 , 33) can 
be calculated limiting attention to the nearest-neighbor 
terms in the summations in Eq. ( p2| ) . We order the three 
principal g factors as 51 < 32 < 33 and parameterize 
them as gj = 2(1 + yj), j = 1, 2, 3. Hence, the tails of the 
distribution correspond \yj\ 3> max(A 2 ,A7y) for at least 
one of the yj. With the definition Q(x) = 1 for x > 0, 
and <d(x) =0 for x < 0, the tails of the distribution are 
found to be 

P(yi <V2< 2/3) = Q(-»a) 



2/3 - yi 

X exp 



(-Trifc)Va 



(yi - 2/2 - 2/3) 2 + 4(yi - y 2 )y 2 

WV2 



(23) 



In the limit 77 A the tail of the distribution factors as 



ut ^ ^ \ J/3 - 2/1 
P{yi < 2/2 < 2/3) oc — exp 

reproducing the result 
3A 2 



2/1 ~ 2/2 

2 



4r] 2 y 2 



P 



47T7/2 



%3)%2 - 2/l) 



(24) 



for the tail of the g tensor distribution obtained in Ref. 
fijlm the limit r/ — > 0. In the opposite limit A 77, Eq. 
(|23|) simplifies to 



f (2/1,2/2,2/3) = 



97 7 2 A 2 9(-2/i)^(2/2)9(2/ 3 ) 
7r (2/i - 2/3) 4 



(25) 



holds if 



Equation (|2J) is valid if |yi|, |y 2 | > A 2 ; Eq. 

M,M > A??. 

In Fig. we have shown the distributions of the princi- 
pal g factors 31, 32, and 33, calculated from the random 
matrix model (Q) using numerical diagonalization. Al- 
though the limits (|5|) and (p5| ) were derived for the tail 
of the (7-tensor distribution only, they can account for 
some qualitative features of the full g-tcnsor distribution 
for weak spin-orbit scattering shown in Fig. |^: when the 
orbital contribution to the g tensor dominates (77 ^> A), 
generically 33 > 2, g 2 w 2, and g% < 2, cf. Eq. (g5h. On 
the other hand, when the Zeeman contribution to the 3 
tensor dominates (77 <C A), one typically has 31 rs g 2 < 2 
and 33 ~ 2, cf. Eq. (§4|). 

We now turn our attention to correlations between g 
tensors of neighboring levels. Such correlations are de- 
scribed by the correlator 



Cij,kl - (gl) 2 



((Gfj,,ijGfj,+l,kl) — (G t j,,ij)(G t j,+l,kl)) ■ 



(26) 



2.0 
1.5 
p (9) 1.0 
0.5 

1.5 
p (9) 1.0 
0.5 
0.0 




9i 

93 




FIG. 

factors. 
\ = g = 



g 

2. Distributions of magnitudes of the principal g 
Upper panel: A = 1.0, g = 0.2; Lower panel: 
0.5. 



Calculating the correlator C to leading order in A < 1, 
we find that the result is dominated by events where the 
levels e M and e^+i are very close. Since this contribution 
is formally divergent, as a result of the presence of the 
energy denominators in the perturbation expression j2^), 
a cut off must be imposed at energy separations £ M +i — 
of order A<5 where the perturbation theory is not valid. 
To treat the contribution from nearby levels e M +i and e M 
correctly, we calculate the contribution from such events 
non-perturbatively. To leading order in A <^ 1, the result 
of such a treatment amounts to the replacement of the 



by l(e 



M-i 



energy denominator e^+i — in Eq. 
e M ) 2 + lApA 2 /^] 1 / 2 , where A is shorthand notation for 
the vector with components , j = 1, 2, 3. We then 

obtain the following result: 



A I 

Cijki = — (SikSji + 6ii5 jk )(r) 2 In \+ — ) 



7T 

3A^ 
IOV 



SijSki- 



(27) 



The correlator between g factors (at a fixed direction 
of the magnetic field) is found from Eq. ( p7f ) setting i — 
j = k = I = B in the direction of magnetic field, 



2A 2 



c = (gi + i9i)/(g 2 r - 1 = — U 2 in a + - 



1 



(28) 



B. Strong spin-orbit scattering 

In the regime of a strong spin-orbit scattering, A > 1, 
the g tensor distribution can be calculated from Eq. ( |2l| ) 
using the known distribution of the eigenvectors of the 
random Hamiltonian (^) at zero magnetic field, see Sec. 
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FIG. 3. g-factor correlation as a function of spin-orbit 
coupling A computed numerically for 200 x 200 GOE-GSE 
crossover matrices. Dashed line shows the result from pertur- 
bation theory Eq. (Eq). 




FIG. 4. Averaged |(j| 2 as a function of spin-orbit strength 
A. The critical value r/o = \/2/3 « 0.81. 



0. We then find that the matrix elements of the 3x3 
matrix G of Eq. (|l|) are Gaussian random numbers, with 
zero mean and with variance 1/A 2 + 2rj 2 . From this we 
coiLclude that the distribution of the principal g factors 



P(9l,92,93) oc 




-3 9l 2 /2<3 2 > 



where 



3 



[9l 



9l 



9l) 



3 

A 2 " 



6r/ 



(29) 



(30) 



Values for r\ for diffusive and ballistic spherical grains are 



given in Eqs. (|Tl| ) and 
tend the result of Ref. 



12ft. Equations (g9j) and (gfj) ex- 
n] to the case rj 7^ 0. Equation 
(30j), which was derived using the random matrix model 
(I|), agrees with the results of Matveev et al, which were 
derived using a comparison of the g factors and the en- 
ergy absorption of a time-dependent magnetic field. 

In Fig. |we show the result of numerical calculations of 
(g 2 ) as a function of the spin-orbit scattering rate A and 
for various values of t). For rj 2 < 2/3, (g 2 ) < 2 for all A, 
while for rj 2 > 2/3, (g 2 ) > 2. The derivatives with A are 
maximal near A = because of the enhanced fluctuations 
due to the orbital part at small A, cf. Eq. (p2|) . 

Correlations between g tensors of neighboring levels 
trivially vanish for large A because, in the GSE, differ- 
ent eigenvectors are statistically uncorrelated. However, 
since the average g tensor also depends on A, it is a more 
meaningful question to study the correlator between g 
tensors, normalized by the average g factor, cf. Eq. (26). 
In the presence of an orbital contribution to the g tensor, 
the average g factors are nonzero for A ^> 1, see Eq. (30), 



so that the vanishing of correlations in the GSE implies 
that they vanish compared to the average as well. With- 
out the orbital contribution, (/-tensor correlations cannot 
be addressed with reference to the eigenvector statistics 
in the GSE, because Q — in the GSE. Instead we need 
the more detailed knowledge of the eigenvector distribu- 
tion for large A, which is summarized in Sec. The 
main result of that section is that the eigenvector distri- 
bution depends on the distribution of certain orthogonal 
invariants p^, j = 1, 2, 3 which are 2x2 matrices in spin 
space, see. Eq. (17). With the help of Eq. (pjl), one eas- 



ily verifies that, in the case rj = 0, the g tensor may be 
expressed in terms of these orthogonal invariants only, 



2tr p'p- 



/-tf-tr f_Lf_L 5 



,7 = 1,2,3, 



(31) 



where the trace is taken in spin space. Since, for A 3> 
1, the orthogonal invariants are all independently 
distributed for different levels, we conclude that g tensors 
of different levels are uncorrelated in the case 77 = as 
well. 

Figure ^ shows the (/-factor correlator (^8|) normalized 
by the average g factor as a function of A. The numerical 
diagonalization confirms our previous conclusions that g 
factor correlations are small for both asymptotic regimes 
A <C 1 and A> 1. Correlations are maximal for interme- 
diate spin-orbit scattering strengths, A ~ 1.5, but never 
amount to more than 10% of the average (g 2 ). 



IV. AVOIDED CROSSING ENERGIES 

Once the Kramers' doublets are split by the mag- 
netic field, half of the levels move upward with slope 



G 



~ (1/2)<7/zb-B, while the other half moves downward with 
the same slope. Hence, a downward moving level e^+i^ 
and the upward moving level e Mi+ meet at magnetic field 
strength 



Br. 



Mb Uu +9n+i)' 



(32) 



In fact, since the matrix element of the coupling Hb to 
the magnetic field between the corresponding eigenstates 
iV-V+i,-) an d lV7t,+) is finite, the two levels do not cross, 
but exhibit an avoided crossing, see Fig. [D. In this sec- 
tion we calculate the minimum distance A between the 
energy levels in the avoided crossing, its dependence on 
the direction B of the magnetic field, and its level-to-level 
fluctuations. 

The avoided crossing energy is well-defined only if the 
magnetic field dependence of the two levels s^+i,-(B) 
and e M) +(-B) is linear, the only exception being the cur- 
vature resulting from their mutual interaction at the 
avoided crossing. For the magnetic field strengths of 
interest, B ~ B c , other sources of level curvature as a 
function of the magnetic field, which arise both from the 
spin and orbital couplings in the Hamiltonian Hb of Eq. 
(H), are small if both A C 1 and r\ <C 1. Hence, for the 
purpose of calculating the avoided crossing energy A it 
is sufficient to consider the perturbative regime of small 
A and small 77. 

Considering the Hamiltonian in the basis of states 
an d 1?/^+), corresponding to the energy levels 
E-ix+i,- and £ M ,+ at zero magnetic field, respectively, 



l i BBg l 



we find that the avoided crossing energy A reads 



A = 2|(^ +1) _|W flc |^ )+ ) 



4k M +i 



\B-M\%,+)\. 



(34) 



Mb (5^+1 +9ti) 
Using first order perturbation theory in A and ?y, we find 



A = A 



B x 



1 



N 



(35) 



plus terms of order Xrj which are not relevant in the 
regime we consider. The components of the vector 
A^ +1 ^ are matrix elements of the spin-orbit matrices 
Aj, j = 1,2,3 of Eq. (^) in the basis that diagonalizes 
the Hamiltonian to zeroth order in A. 

In order to find the distribution of the avoided crossing 
energy A, we write 



A = An sin t 



(36) 



where < 9 < tt is the angle between the direction B of 
the applied magnetic field and the vector A^ +1 ^ . Using 




2 

A/<A> 

FIG. 5. Main panel: Distribution of the avoided crossing 
energy A. Solid line is the perturbative result (p9|) ; the data 
points are from numerical evaluation of Eq. (BJ) using the 
numerical diagonalization of the random matrix model (^) 
with r) — and A = 0.2 (crosses), and A = 0.4 (circles). 
Inset: Comparison of perturbation theory (solid curve) and 
numerical results (data points) for the average (A). 

the known distribution (0) of the spin-orbit coupling ma- 
trices Aj (j = 1, 2, 3), one finds that the three elements of 
^p+i-p eacn have a Gaussian distribution with zero mean 
and with variance N5 2 /tt 2 . Hence, we conclude that the 
vector A^ +1 '^ is randomly oriented in space, so that 



P(6) = - sin9, 



and that 



(Ao^)V2^ 
(XSf 



■ exp 



vrAo 
XS 



(37) 



(38) 



Equations (36)-(|38|) not only give the full distribution of 



the avoided crossing energy A, but also the dependence 
of A on the direction B of the magnetic field. Equations 
(BQ)-(p8h can be combined to give 



P(A) 



7T 2 A 

Jxsy- 



■ exp 



1 /ttA 

2 (U 



(39) 



The latter result is relevant for comparison with exper- 
iments where the direction of the magnetic field cannot 
be variedBQ 

Figure |B| shows the distribution (|3^) , together with re- 
sults from a numerical calculation of the distribution of 
Eq. ([34]) using the random matrix model (^) for r\ = 
and two different values of A. We see that the agreement 
between the numerical diagonalization of the random ma- 
trix model and the distribution ( |39"| ) calculated using first 
order perturbation theory in A remains good up to A ~ 1 . 
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[We should note, however, that the approximations lead- 
ing to an avoided crossing energy that is dominated by 
matrix elements involving two neighboring levels only, is 
valid for A <C 1 and r\ <C 1 only, see the discussion pre- 
ceding Eq. (|34|).] Although there are corrections to P(A) 
to second order in A, the first nonzero corrections to the 
average (A) appear to third order in A only. 



V. DISCUSSION AND CONCLUSIONS 

In this paper, we have presented a random matrix the- 
ory for the distributions of g tensors and avoided crossing 
energies in small metal grains with spin-orbit scattering. 
Our theory includes both the spin and the orbital effects 
of the magnetic field. 

For large spin-orbit scattering, the main effect of the 
orbital contribution is to increase the typical size of the g 
tensor; the fluctuations (normalized by the average) and 
the relative magnitudes of the three principal g values axa 
the same with and without a large orbital contribution.!!!! 
For weak spin-orbit scattering, the presence of an orbital 
contribution to the g tensor not only increases the av- 
erage of the g-tensor distribution, it also changes the 
relative magnitudes of the principal g values. Without 
orbital contribution, two principal g values are approxi- 
mately equal and smaller than two, while the third prin- 
cipal g value is close to 2. If the orbital contribution is 
large, all three principal g values are different and, on 
average, symmetrically positioned around two. 

Petta and Ralph have measured distributions of g fac- 
tors (i.e., the square root of the Q zz element of the g- 
tensor) for small particles of different metals and found 
that distributions, if normalized to the average, were in 
very good agreement with the random matrix theory of 
Ref. fLU The average of the distribution, however, was 
up to a factor 10 smaller than the theoretical predictioe, 
( ^0| ) with a reasonable estimate for the parameter rjS3 
A similar discrepancy between a experimental and the- 
oretical estimates was reported in a different context by 
Marcus et a/.Ej for the magnetic field scale for fluctu- 
ations of Coulomb blockade heights in two-dimensional 
/j,m-size GaAs/GaAlAs quantum dots (see also Ref. 21). 
Although the experimental system studied in Refs. |20|j2f| 
is quite different from that of Petta and Ralph, the ran- 
dom matrix theories describing the magnetic field depen- 
dence of Coulomb blockade peak heights and the orbital 
contributions to g factors are the same. At present, we 
do not know of a solution to either puzzle. 

One complication in the search for an orbital contribu- 
tion to the g factors measured in Ref. ^ is that the main 
effect of the orbital contribution is to change the average 
of the g-factor distribution only. Since, for strong spin- 
orbit scattering, the average g factor depends on both the 
dimensionless spin-orbit coupling A and the dimension- 
less orbital contribution 77, cf. Eq. (|30|), it is impossible 
to characterize what fraction of a measured g factor is 



the result of a state's orbital magnetic moment. The re- 
cent development j-pf experimental methods to measure 
the entire g tensora opens new avenues to investigate the 
orbital contribution. For weak spin-orbit scattering, the 
<?-tensor distribution depends on the two parameters A 
and r\ in a nontrivial way; even a weak orbital contri- 
bution leads to g tensors with, at least, one principal 
g value larger than two, see, e.g., Fig. Hence, mea- 
surement the full g tensors for metal grains with weak 
spin-orbit scattering, such as large Al grains, eventually 
doped with a small concentration of Au,H will allow the 
independent determination of the orbital contribution. 
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